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PACS 64 . 60 . Ht - Dynamic critical phenomena 

PACS 64 . 60 . Kw - Multicritical points 

PACS 05. lO.Cc ~ Renormalization group methods 

Abstract. - A complete two loop renormalization group calculation of the multicritical dynamics 
at a tetracritical or bicritical point in three-dimensional anisotropic antiferromagnets in an ex- 
ternal magnetic field is performed. Although strong scaling for the two order parameters (OPs) 
perpendicular and parallel to the field is restored as found earlier, in the experimentally accessible 
region the effective dynamical exponents for the relaxation of the OPs remain different since their 
equal asymptotic values are not reached. 



' ^ [ Systems with more than one order parameter (OP) ex- 
^ hibit a rich variety of phases separated by transition lines 
^ which might meet in multicritical points. The interaction 
' b light favor simultaneously ordering of two OPs. Such a 
doubled ordered phase is known as supcrsolid phase [T] 
^ and is under investigation since its possible observance in 
^^'*He [2]. The two OPs might describe the order in phys- 
ically different phases: E.g. from the normal fluid to the 
^uperfluid or to the normal solid phase. If both orders 
[ appears one might have transitions from the superfluid to 
^^Othe supersolid phase and from the normal solid phase to 
^^the supersolid phase. Another example would be a system 
With transitions to a superconducting and a magnetically 
• 'Ordered phase and a phase where both orderings appear. 
. ^ In purely magnetic systems the phases are characterized 
by different orderings in spin space. There is a correspon- 
^—^ dence between the quantum liquid system and magnetic 
systems where the supersolid phase corresponds to the bi- 
conical phase [3] . The existence of a biconical phase leads 
to the occurrence of a tetracritical point where four second 
order phase transition lines meet and which belongs to a 
new universality class [3]. 

In the case of the three-component (n = 3) three- 
dimensional {d = 3) anisotropic antiferromagnets in an 
external magnetic field in z direction the disordered (para- 
magnetic) phase is separated from the ordered phases by 
two second order phase transition lines: (i) one to the spin 
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flop phase (ordering in the spin space perpendicular to the 
external magnetic field) and (ii) one to the antiferromag- 
netic phase (ordering parallel to the external field). The 
point where these two lines meet is a multicritical point 
which turned out to be either tetracritical or bicritical de- 
pending on whether the ordered phases are separated by 
an intermediate biconical phase or not. The static phase 
transitions on each of the phase transition lines belong for 
(i) to an XY-model with n = 2 and for (ii) to an Ising 
model with n = 1 [J] . Whether a bicritical or tetracritical 
point is realized depends on the specific fourth order cou- 
plings HUni . Several antiferromagnets have been suggested 
for observing the theoretically proposed phase diagrams, 
a review on the experimental situation can be found in [S]- 

Concerning the dynamical universality classes they 
might be different for the systems mentioned above de- 
pending on the possible reversible and non reversible terms 
in the equations of motion and the conservation properties. 
In the magnetic system considered here the transition (i) 
belongs to the class described by model F and (ii) belongs 
to the model C class (for the definitions of the models 
see [9]). At the multicritical point the critical behavior is 
described by a new universality class both in statics and 
dynamics characterized by the biconical fixed point [S]. 
The advantageous feature of these systems is that all the 
different OPs characterizing the ordered phase are physi- 
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cally accessible. This is most important for the dynamical 
behavior since the only other example belonging to model 
F is the superfluid transition in "^He where the OP is not 
directly measurable. Here the OPs are the components 
of the staggered magnetization. Their correlations (static 
and dynamical) are experimentally accessible by neutron 
scattering. Realistic models might be more complicated 
(see e.g. [TU]) but the behavior near the multicritical point 
is well described by the rcnormalization group (RG) the- 
ory. 

The dynamical model we analyze goes beyond the pure 
relaxational dynamics [TT] and has been considered by 
means of the field theoretical RG approach in [T^HTl] re- 
placing earlier mode coupling theories [TS]. It was argued 
that due to nonanalytic terms in e = 4 — d a dynamical 
fixed point (FP) in two loop order (which was calculated 
only partly) qualitatively different from the one loop FP 
is found. In one loop order the relaxation times of the 
components of the staggered magnetization parallel and 
perpendicular to the external magnetic field scale differ- 
ently whereas in two loop order they would scale similarly 
if the new FP would be stable. In addition it turned out 
that the FP value of the timescale ratio of the two OPs 
cannot be found by e expansion and might be very small 
at d = 3, namely of 0(10"^^). A basic assumption of 
the above analysis was that within statics the Heisenberg 
FP is stable. However it turned out in two loop statics 
using resummation techniques that in d = 3 the Heisen- 
berg FP interchanges its stability with the biconical FP 
[5]. Here we calculate the complete functions in two loop 
order which allows us to consider the non-asymptotic be- 
havior near the multicritical point. 

The non-conserved OP in an isotropic antiferromagnet 
is given by the three-component vector <f)Q of the stag- 
gered magnetization, which is the difference of two sublat- 
tice magnetizations. In an external magnetic field applied 
to the anisotropic antiferromagnet the OP splits into two 
OPs, (/)_Lo = ((/)§, (jy^) perpendicular to the field, and 
'/'II = 4>o parallel to the external field. In addition to the 
two OPs the z-component of the magnetization, which is 
the sum of the two sublatticc magnetizations, has to be 
considered as conserved secondary density mo. The static 
critical behavior of the system is described by the func- 
tional 
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{r} and field h [SJ[TT]. One may switch from the de- 
scription of real OP components (/jj^q to a complex OP, 
a macroscopic wave function, as it appears in a superfluid 
or superconductor by defining tJjq — (j)Q ~ i(f>Q . Apart from 
demonstrating that all these systems belong to the same 
the static universality class it also is of practical advantage 
in the dynamic calculation. 

The critical dynamics of relaxing OPs coupled to a dif- 
fusing secondary density is governed by the following equa- 
tions of motion [T2] (there the complex OP tpo was used): 
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with the Levi-Civita symbol e''^'^. Here a,/3 = x,y and 
the sum over repeated indices is implied. Combining 
the kinetic coefficients of the OP to a complex quantity, 
T±_ = T'j_ + ir" , the imaginary part constitutes a preces- 
sion term created by the rcnormalization procedure even 
if it is absent in the background. The kinetic coefficient A 
and the mode coupling g are real. The stochastic forces 
dd,. , ddi. and 6m fulfill Einstein relations 
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The reversible terms in these dynamic equations for the 
OP components and the conserved density have been de- 
rived by using generalized Poisson brackets for the spin 
components defining the staggered magnetization and the 
z— component of the magnetization (for more details see 
[9] and the literature cited there). An exception consti- 
tutes the term with F" which appears due to the rcnor- 
malization procedure and the nonzero asymmetric cou- 
pling tI. Similar dynamic equations may also appear in 
the other systems where biconical phases are observed. 
The superfluid transition is described by model F and also 
for the superconducting transition this universality class 
has been suggested [TOlfTT] although no explicit derivation 
based on the methods used here has been performed. 

Applying the rcnormalization procedure using minimal 
subtraction scheme jl8j we find the fiow equations for the 
time scale ratios of the renormalized kinetic coefficients 
and the mode coupling between the perpendicular OP 
components and the magnetization. We define time scale 
ratios by the ratios of the kinetic coefficients of the OPs 



with familiar notations for bare couplings {u,j}, masses and the secondary density w± 
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as the ratios between the relaxation rates of the two OPs 
V = fT- = = ^ = and the mode couphng pa- 



rameters f± = g/ ^r'j_A or F = g/X. For these dynamic 
parameters we obtain the flow equations 
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where I is the RG flow parameter and the ^A-function is 
obtained by the rcnormalization procedure as 
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The function Q = Q{j±,w±, F) contains all higher order 
contributions beginning with two loop order and is identi- 
cal to the corresponding function in model F (see (A. 28) 
and (A. 29) in [9]). We obtain the C- function for the per- 
pendicular kinetic coefhcient F^ as 
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where we have introduced the coupling D± — w±"f± — 
iF. The functions A_i = Aj_{^^,Tj^,w_i, F), B± = 
B±{'j±,T±,w±, F) are identical to cqs. (A.25), (A. 26) in 
[9]. X± is defined as 
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Cp'^^ ({u}, w^, w) is the ^-function of the perpendicular re- 
laxation Fj^ in the biconical model A, but now with a 
complex kinetic coefficient Fj^ 

C(lHM,.x,.) = f (21n^ 
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The dynamic ^-function of the parallel relaxation kinetic 
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coefficient F|| is obtained as 
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Table 1: Two loop FP values of the mode coupling f±, the 
ratios q — uiy/wYi s = w'^/w'j_ and the dynamic exponents 
in the subspace tun = 0, w± = with finite value of v = 
q/{l + is) for the static biconical B and Heisenberg FPs. 
For comparison we add the FP values for the exponents that 
govern critical dynamics at magnetic fields below and above 
the multicritical point. These are described by model C at 
n = 1 and model F at n = 2. 
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~(C) (A*) 

Cr,, (""IhTih^;!!) = Cr(tt|| , 7|h T]! , u'n) - Cr (""Ihril)' '^'here 
the functions on the right hand side are defined by (A. 8) 
and (A. 9) for n = 1 in [5]. The functions Ti and T2 are 
defined as 
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and Cl^'' ({u},v±,v) is the ^-function of the kinetic coef- 
ficient of the parallel relaxation in the biconical model A. 
With a complex F^ it reads 
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In order to find the FP values of the time-scale ratios 
and the mode coupling the right hand sides of eqs. ([8|)-(l9]) 
have to be zero. If the FP value of the mode coupling f± 
would be zero one obtains the FP values of the time ratios 
of model C discussed in [Tn| . However this FP is unstable. 
If the FP value of f± is nonzero then due to the logarithmic 
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Fig. 1: Effective dynamic exponents in the background using 
the flow equations ((H]),® in the complete dynamical parame- 
ter space. The static values are taken for the Heisenberg FP 
(dashed curves) and for the biconical FP (solid curves). 

terms in the ^-functions both OP have to have the same 
time scales i.e. a finite nonzero FP value v* . This is only 
possible either for nonzero finite FP values of and w\\ 
or when both of these FP values are zero. No finite FP 
values for wj_ and have been found. In the other case 
the approach to zero of both time scales has to be the same. 
Therefore the approach to the multicritical dynamic FP is 
described by the flow in the limit w± 0, w\\ and 
V finite (asymptotic subspace). The flow in the complete 
dynamic parameter space and in this asymptotic subspace 
will be discussed afterwards. 

The ^-function for the perpendicular OP relaxation 
might be complex, Cr^ = Cf + *Cf ■ order to ob- 
tain the usual asymptotic power laws for the relaxation 
coefficients Fn and F^ the FP value of the imaginary part 
C,"* has to be zero. In consequence the asymptotic flow of 
the real and imaginary parts of v is governed by the same 
exponent - C?,, ■ 

If the FP value of the mode coupling / is different from 
zero and finite one has from eq. (O e + Cf*^ -I- c}^^* = 
and the relation [M] between dynamical and static critical 
exponents z± + Zm = 2^ (here the z exponents govern the 
corresponding scaling times and and v are the crossover 
and correlation length exponents). The dynamical expo- 
nents are defined as Zq = 2 + Co with o =±, ||, m. Because 
V* is finite and nonzero z±_ = = ^op- This means 
that strong scaling with respect to the OPs, the compo- 
nents of the staggered magnetizations, but weak scaling 
with respect to the conserved density, the magnetization 
TO, holds since Zm ^ ^^op- 

The two loop order values of the dynamic exponents 
together with the FP values of the time scales and the 
mode coupling are presented in table \T\ For the model of 
the three-dimensional anisotropic antiferromagnet under 



Fig. 2: Effective dynamic exponents in the asymptotic sub- 
space ttJii — wj_ — and v = 'w^^/w± ^ and finite. Dashed 
and solid curves as in Fig. [1] 

consideration, the biconical FP B (u\ 7^ U| ^ u*) has 
been shown to be stable. It governs the static critical 
behaviour in the complete space of couplings (see e.g. [S|). 
Substituting their two-loop values obtained in [S] within 
generalized Pade-Borel resummation technique |22| into 
the fiow equations (IHl)-(|ni), we get two dynamical FPs, 
their coordinates are given in the first and third row of 
table [T] Although two different dynamical FPs arc found 
(with zero and nonzero s*) this difference does not lead to 
a change in the corresponding FP values of the dynamical 
exponents. This is because both FPs have extremely small 
but different g*. 

For comparison we have included, besides the biconi- 
cal FP B describing tetracritical behavior, the isotropic 
Heisenberg FP H (u\ = u| = u* ) describing bicritical 
behavior. This FP is only reached in the subspace of the 
static couplings that lie in its attraction region (see fig. 3 
in [S]). Again, as in the case of the biconical FP B^ we 
obtain two dynamical FPs with coordinates given in the 
second and fourth line of table [1] The numerical values of 
the dynamical exponents are practically equal in the dif- 
ferent dynamical FPs Wc further quote in table [T] the 
dynamical critical exponents on the two phase transition 
lines below and above the multicritical point, which are 
given by model C and model F respectively. 

The FP value of v is extremely small and therefore in the 
physical accessible region one cannot prove strong scaling 
for the OP components. Indeed in the non-asymptotic 
region the dynamic parameters are described by the flow 
equations ^ , ([9]) and from these dependencies the effective 
dynamic exponents can be calculated. The result is shown 
in fig. [T] The static parameters have been set already to 
their FP values and therefore the starting values of the 
effective exponents are different from z = 2. It turns out 
that the pref actor of the Inw-terms in eqs. ([TT|) and ([T^. 
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which drive the flow of the dynamic parameters into the 
asymptotic subspace is reduced and the flow is almost like 
in one loop order. Therefore weak scaling with z^^^^ ^ 2.04 

and z^l^ ^ 1.6 is observedQ 

The approach of the effective dynamical exponents in 
the asymptotic subspace = W|| = and v finite to their 
biconical FP values is shown in fig. [H The background 
behavior is dominated by a behavior corresponding for the 
perpendicular components by model F and for the parallel 
components by model A with a finite value Re(v) whereas 
Im(w) is almost zero. Therefore for the biconical case even 
for flow parameter values InZ ~ —10^ the two effective 
exponents do not reach their asymptotics: z^^^ < zqp and 
z^-^ > Zra- This is different for the Heisenberg case where 
the FP values of the dynamical exponents are reached (see 
dashed curves in fig. . 

Although from our calculation we conclude that the 
asymptotics (strong scaling) would be unobservable effec- 
tive exponents as described in fig. [1] are observable. The 
complete two loop calculation allowed us to calculate not 
only the values of the dynamic FP but also the effective ex- 
ponents which are the quantities governing the behavior of 
the transport coefficients, i.e. the relaxation and diffusion 
coefficient of the staggered magnetization and magnetiza- 
tion respectively. It is well known that near a dynamical 
stability boundary separating a strong scaling FP with a 
finite timescale ratio from a weak scaling FP with a van- 
ishing time scale ratio also small dynamic transient expo- 
nents appear and effective critical behavior is observed. 
The case where the OPs have the component values n = 2 
and ?i = 1 is located near the stability boundary between 
the biconical and decoupling FP as has been demonstrated 
in [5]. For the decoupling FP the time scales of the two 
OPs scale differently and weak scaling is expected. The 
effective values of the dynamical critical exponents (start- 
ing from different initial conditions) are driven to almost 
stationary values. These might be measured in neutron 
scattering experiments. 

A natural question concerns the reliability of numerical 
predictions for the observables obtained in our study. In 
particular, how will an increase of the order of the per- 
turbation theory influence the numerical estimates. An 
estimate can be given by comparing results obtained in 
different perturbation theory (loop) orders. This can be 
done for the static part of the RG functions, which are cur- 
rently known with a record five loop accuracy [23]. As it 
was demonstrated in Ref. [S], the two loop approximation 
we consider here refined by the resummation is enough 
to catch the main features of the static phase transition 
(FP stability and respective universality classes) as well 
as to give reliable estimates for the observable quantities 
that govern the phase transition (leading exponents and 



^Even for flow parameters In/ —10® no visible changes in the 
values of the different z'^^ ^ occur although the dynamic parame- 
ters change. Thus the asymptotic subspace is not reached for these 
extremely small values of I. 



corrections to scaling) . It is well known, that the dynamic 
RG calculations are technically much more complicated as 
static ones. In particular, no higher orders of the pertur- 
bation theory are known for the model we consider here. 
However, as known from the previous experience in the 
RG description of dynamical criticality [5] we expect also 
in this case that the two loop calculation captures the es- 
sential dynamical properties, namely to be effectively in a 
weak scaling situation. 

The comparison between experiment and theory for 
multicritical behavior is much less developed in dynam- 
ics than in statics (see e.g. the situation at the tricriti- 
cal point |24|). But exploring even the critical dynamics 
along the two transition lines is of interest since the OP 
and the condserved density are experimentally accessible. 
Also computer simulation might be considered for a com- 
parison |25] and explicit theoretical results arc worthwhile 
for the interpretation of the numerical results. Not only 
exponents are necessary for a careful interpretation but 
also the calculation of the dynamic structure factors arc of 
interest. This also concerns the neutron scattering experi- 
ments. The effective values of the timescale ratios, known 
from flow equations (like (|H1)), enter the shape functions 
and may change their shape from a Lorentzian consider- 
ably. 

^ ^ ^ 
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